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In this Letter we study discrete symmetries of mean field manifolds of condensates of F — 2 
cold atoms, and various unconventional quantum vortices. Discrete quaternion symmetries result 
in two species of spin defects that can only appear in integer vortices while cyclic symmetries are 
found to result in a phase shift of 2n/3 (or 4ir/3) and therefore 1/3- (or 2/3-) quantum vortices in 
condensates. We also briefly discuss 1/3-quantum vortices in condensates of trimers. 



One of the most striking features of a superfluid is the 
existence of quantized vortices. This is a consequence of 
the requirement that the quantum mechanical wavefunc- 
tion of any physical state be single- valued. In a standard 
single component bosonic condensate, this requirement 
results in the quantization of circulation: the line-integral 
of the superfluid velocity along a path linking a vortex 
takes discrete values: § dr ■ w s — n ■ with n an in- 
teger. Here, is the fundamental quantum of circula- 
tion. The quantization of circulation has been observed 
in superfluid liquid Helium by Vinen [l[ , and recently in 
Bose-Einstein condensates of cold atoms Q ; the presence 
of quantized vortices usually results in discrete values in 
the energy splitting of collective modes that have been 
studied in experiments. 

The existence of quantized vortices in a condensate is 
also tied to the topology of the vacuum manifold. In 
a single component condensate, the vacuum manifold is 
the set of all distinct phases that the condensate can 
have. This set is identical to the unit circle S 1 . As one 
moves in a closed path around a quantized vortex,the 
phase must traverse S 1 an integer number of times [H, 0, 
[H, For condensates with more than one component, 
vortices can have circulation quantized as fractions of the 
basic unit ; this is achieved when a change of phase is 
combined with a rotation of spin or orbital orientation. 
A most recent example are the half-quantum vortices in 
condensates of sodium atoms 0; half vortices were also 
discussed in the context of superfluid 3 iJe[3|. 

Cold atoms in optical traps and lattices are a promis- 
ing venue where exotic condensates and vortices could be 
realized in nature. Spin correlated states of spin F = 1 
cold atoms have alr eady been explored experimentally Q 
and theoretically [l,[l(J, 11|. More recently, the scattering 
lengths of spin F = 2 cold atoms have also been inves- 
tigated and various possible condensates have been 



pointed out ID, 14\. In this letter we shall examine the 
possible appearance of unconventional vortices in the lat- 
ter condensates. In a previous study of their possible in- 
sulating phases [l5[ , we found a convenient parametriza- 
tion of the F = 2 wave-function as a tensor Xa/3 where 
a,fl take the values x,y,z. Being the wave-function of 
an atom, the entries of this matrix are complex func- 
tions. To be a state with F = 2, it must be symmetric, 



Xa/3 = X/3a, and traceless, tr% = Xxx + X yy + Xzz = 0. 
This leaves five independent complex entries. These en- 
tries are linear combinations of the five components of 
the atom wave-function ip mF with tuf = 2,1,0,-1,-2 
the z-component of the spin. In Ref . [lj| , we related the 
Landau-Ginzburg free energy functional of the conden- 
sate to recently measured scattering lengths for atoms. 
By including a derivative term, in the long-wavelength 
limit we have 



drp [—trV X *(r) ■ V*(r) + MA( X (r),x*(r))], 

TO* 



A( X ,xl = ±b L Tr[ X ,xT + 2c L tr X 2 Tr X 



(1) 



to* is the effective mass of atoms in optical lattices, M is 
the average number of atoms per site and po the average 
particle density. The coefficients br, and cr, can be esti- 
mated from the two-atom interaction strengths in chan- 
nels with various spins [l^; br, or cl is proportional to 
b — (a^ — a-i)jl or c = (a® — a±)jh — 2(a2 — a,±)/7 and ao,2,4 
are two-body s-wave scattering lengths in F = 0,2,4 
channels respectively. For ground states we consider uni- 
form condensates. The x _ dependence of the energy of 
condensates is identical to that of Mott states obtained 



earlier in Ref. 15[; minimization of this energy subject to 
a constraint of trx*x =1/2 was previously carried out. 

When all scattering lengths are equal, A(x,X*) = 
and the Hamiltonian is SU (5) symmetric. The mean field 
manifold is U (1) x SU(5). SU(5) is simply connected and 
its fundamental group is trivial. Only integer vortices are 
allowed in this case. Restoring the spin-dependent scat- 
tering turns on the potential and, as shown below, results 
in multiply connected mean field manifolds and exotic 
vortices. Here we investigate various discrete symmetries 
of the manifolds for cyclic and nematic phases. 

When bh > 0, and cl > 0, the energy is minimized by 
a ground state x satisfying tr% 2 = = [x, X*]- The real 
and imaginary parts of x commute with each other and 
we can find a coordinate frame where x is diagonal. A 
solution which obeys tix = as well as t?x*X = \ i s 
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The full set of degenerate solutions can be obtained by 
applying symmetry transformations: SO (3) rotations, 
implemented by 3 x 3 orthogonal matrices 1Z and mul- 
tiplication by a £7(1) phase, e^, to the solution x+- An 
arbitrary solution \ can be written as e^TZx+TZ -1 ■ 

The set of all solutions is equal to the group of all 
such symmetry transformations factored by the subgroup 
which leaves x+ invariant. As we shall explain below, this 
subgroup is the tetrahedral group, T. The vacuum man- 
ifold is therefore the set of right cosets M. — s °( 3 )^- u ( 1 ) 
[ill]. The tetrahedral group, T, is the subgroup of 
50(3) x £7(1) which leaves x+ m e q- © invariant. This 
group contains the identity, and three diagonal SO (3) 
matrices (called the Klein 4-group) 
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The remaining elements of T implement cyclic permu- 
tations of the diagonal elements of x+ as weu as multi- 
plication by phases. The the elements of SO (3) which 
permute the diagonal elements of x+ are 
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Indeed, x+ = ^Cx+C -1 , where u> was introduced in 
Eq. ([2]). The elements 1,C,C 2 form a cyclic subgroup 
of SO (3), C3. We emphasize that the transformation of 
X+ involves a U(l) phase of ^ as well as a cyclic rota- 



tion represented by C. We will denote this combination 
of transformations as C. The full tetrahedral group T 
has the 12 elements T — {l,I x ,X y) X Z) C,l x C,l y C,T z C, 

C J-^xC l^yC ^X Z C }. 

The first step in classifying vortices is to identify the 
homotopy group IIi(A / J)[l7]. Since 50(3) is not simply 
connected, to understand the topology of the manifold 
M we shall lift 50(3) to 5£7(2). An element of SU(2) is 
a 2 x 2 unitary matrix. The set of 2 x 2 unitary matrices 
is characterized by the four real Euler parameters (eo, e) 
through the relation Q = eo + ie ■ B\ a are the usual 
Pauli matrices. Q is unitary when eg + e • e = 1; the 
Euler parameters live on the unit three-sphere 5 3 , which 
is simply connected. The matrix Q is the well-known 
quaternion representation of a rotation [l8| . The 3x3 
rotation matrix 1Z € 50(3) which corresponds to Q e 
SU(2) is IZap = ^Tr (caQa^Qt). This is a two-to one 
mapping since both Q and — Q are mapped to the same 
TZ. The inverse, 1Z -> (Q, Q) is the 'lift' of 1Z to SU(2). 

We also need to lift the tetrahedral group T to T*, the 
binary tetrahedral group which is a subgroup of SU (2) 



and contains 24 elements. For this, we must identify the 
pair of SU (2) matrices which correspond to each of the 12 
elements of T. The elements of T are the simple rotations 
(examples are shown in Fig© that lift as: 1 — > ±1, 



-x,y,z 



±i<7x,y,z, C -* ±<J = ±|(1 



ia z ). 



We further denote a rotation a combined with a phase 
shift of 27r/3 as a. The group elements of T* which are 
lifts of T are therefore, 

±ia x , ±i(T y , ±i<r z , ±1 (5) 
±ia x a, ±i<r y a, ±ier 2 (T, ±cr (6) 
iiiTaCr 2 , ±ia y a 2 , ±ia z a 2 , ±cr 2 . (7) 

Eq.© represents the quaternion subgroup Q. 

The vacuum manifold derived above is isomorphic to 



M 



SU{2) x U(l) 



(8) 



To find the homotopy classes of closed paths in M., we 
must first understand how they are related to paths in the 
product of space of SU(2) x U(l). As we have discussed 
above, SU(2) is equivalent to a three-sphere 5 3 , and U(l) 
to a circle 5 1 , so that we could equivalently think of paths 
on the product space 5 3 x 5 1 . Each point of 5 3 x 5 1 
corresponds to an element in the group SU(2) x £7(1). 
Factoring by T* is simply making an identification on 
5 3 x 5 1 , each point is to be identified with 23 other points 
which are obtained from the corresponding group clement 
by operation of the non-trivial elements of T* . 

Consider a path on 5 3 x 5 1 . In order to be a closed 
path on M it must either end at the point on 5 3 x 5 1 at 
which it began, or it must end at one of 23 points that 
are identified with it. As shown in Fig. [l] there are three 
types of paths (which all begin at the identity): 

Type A: Paths that traverse 5 1 n (n is an integer) 
times and end with one of the eight elements in Eq.©; 

Type B: Paths that traverse 5 1 n + | times and end 
with one of eight elements in Eq.©; 

Type C: Paths that traverse 5 1 n + | times and end 
with one of eight elements in Eq.©. 

Thus each path is characterized with two variables: 
the winding number around 5 1 and an element in 
T*. Hi (M) = {(n,l),(n,-l),(n,ia a ),(n,-ia a ), 
(n + 4)5-)) (n -I- |, — cr), (n + ^,ia a a), (n + |, —ia a a), 
(n + |,CT 2 ),(n + |,-CT 2 ),(n+ j,ia a a 2 ),(n+ ^,-ia a a)} 
This set forms an infinite order non-abelian dis- 
crete group with composition law (x,g±) ■ (^,32) = 
{x + y,gig 2 )- 

So far, we have computed the homotopy group IIi(A^) 
for loops with a base-point. However, distinct vortices 
are identified with free homotopy classes, which coincide 
with conjugacy classes of the based homotopy group 4, 
5]. For an element gi of a group, the conjugacy class 
is the set of all distinct elements in the set g a 99a l as 
g a sweeps over the group. In our case, the classes are 
1) {(X 1 )}; 2 ) 3 ) {(n,ia a ),(n,-ia a )}, a = 
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FIG. 1: (color online) Schematic of the mean field manifold 
of cyclic condensates. The unit circle is oriented horizontally 
with phase angles £ (defined after Eq(2} around it shown ex- 
plicitly. Each cross section at a given angle stands for a three- 
sphere. Examples of 24 four identified points in the product 
space S 3 x S 1 are shown explicitly as stars, with eight points 
at £ = or 2-7T (or elements in Eq(5j , eight points at £ = 2n/3 
(or elements in Eq(6]), and the eight points at £ = 4ir/3 ( or 
elements in Eq[7J). The red path ending at phase 2n is for an 
integer vortex with a spin defect inserted; the red dashed one 
is for a plain integer vortex without spin defects. The short 
(long) blue path ending at 2n/3 (47r/3) is for a 1/3 (2/3)- 
quantum vortex. The green path at phase £, — corresponds 
to a spin defect and the winding along S 1 is trivial. 



x,y,z; 4) {(»+|,cr), -ia a a)}; 5) 



ki<r a a)}; 6) {{n + §,c^j,(n + ^-ia a a 2 )}; 7) {(n + 
|, — a 2 ), in + |, ia a a 2 y\. Altogether there are only seven 
distinct classes and thus seven types of linear defects. 

If n — 0, a Type- A path corresponds to a pure spin 
vortex. Such a path might begin at the origin (iden- 
tity) and end at one of ±io~ x , Ho~ y , Ha z or —1. These 
correspond to, respectively, 180° rotations around one of 
the three axes x, y, z or a rotation around the z-axis by 
360° (see (b-e) in Fig. @). However, paths which end 
at Hu x , ±i<7 y , ±io~ z are freely nomotopic to each other 
and all represent the same spin defect. So there are two 
differ types of spin vortices (specified by classes): (0, —1) 
and (0,±ia a ) 7 a = x, y or z. 

On the other hand when n is a nonzero integer, a type- 
A path corresponds to an integer vortex. Three conju- 
gacy classes involving integer n represent three distinct 
species of integer vortices: a conventional integer vor- 
tex with no spin structure, an integer vortex with a spin 
defect of (0,-1) inserted and an integer vortex with a 
spin defect of (0, ±ia a ) inserted. The last two types are 
unique to condensates of F = 2 atoms. 

Along Type-B or Type-C paths, a phase shift of 27r/3 
or 47r/3 is generated under the action of the element a c 
or a 2 . In these cases, the spin rotation leads to a 1/3- 
quantum vortex or a 2/3-quantum vortex. Furthermore, 
there are only two topologically distinct species of each 
of 1/3- or 2/3-vortices. The two distinct 1/3-quantum 
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FIG. 2: Spin rotations in integer vortices, 1/3-quantum vor- 
tices and 2/3-quantum vortices, (a) the orientation of three 
basis vectors xyz at azimuthal angle <f> = 0. (b-e) the orienta- 
tion of the triad xyz at angle cj> — 2tt in four configurations of 
integer vortices: (0,ia x ) (b), (0,ia y ) (c), (0, ia z ) (d), (0,-1) 
(e); the arrow in (e) represents a 360° rotation around the z 
axis. (f,g) the orientation of the triad xyz at angle <f> — 2n 
in two 1/3 quantum vortices: (1/3, ia z a) (f), (1/3, ia y a) (g); 
(h,i) the orientation of the triad xyz at angle cj> = 2-7T in two 
2/3 quantum vortices: (2/3, i<7 z <r 2 ) (h), (2/3, ia x a 2 ) (i). 



vortices are represented by either paths that end with 
elements a, —io~ a a or paths that end with —a, icr a a. 
The 1/3 and 2/3-quantum vortices necessarily contain 
topologically non-trivial spin configurations which differ 
in structure from those in integer vortices. 

We now are going to examine spatially slowly varying 
matrix x(r) and focus on spin vortices and 1/3-quantum 
vortices. All linear defects are oriented along the z- 
direction; the center of a defect is at the origin of cylin- 
drical coordinates (p, 4>, z). The superfluid velocity in the 
condensate is v s = 21m trx*Vx^-/™. Far away from the 
center, the wavefunction of a pure spin vortex is speci- 
fied by the rotation, \(p = +00, 0) = IZ^x+Tt 7 '(</>) . To 
minimize the energy in Eq.([T]), we find that the rotation 
matrix satisfies the equation 



[X+, -rrr] = 0, where A = Tl T — . 



(9) 



A solution can be expressed as TZ(<j)) — exp(in • L/(0)), 
where f(4>) is a linear function of </>. Alternatively 

Tla0(<t>) = <5 Q /? cos /(<?!>) + sin/(0)e Q/ 3 7 n 7 

+n a n /J (l-coB/($). (10) 

Here L is a matrix vector, = — ie Q /3^, e Q /3 7 is the 
antisymmetric tensor, n is a unit vector with three com- 



4 



ponents n a , a = x,y, z. For a configuration that is spec- 
ified by the boundary condition: 1Z(0) = 1,1Z(2tt) = l z , 
we find the following solution 

/(^) = |,n = e z . (ff) 

The superfluid velocity is zero. When going around this 
vortex, the local triad of three orthogonal basis vectors 
xyz makes a f80 rotation around z-axis (see Fig.([2}l) ). 

In a I/3-quantum vortex, the boundary condition is 
n{<t> = 0) = 1, n(<t> = 2tt) = 1 Z C; at any 0, x(<t>) = 
e^/ 3 ^(0)x+^ T (0). The solution which satisfies Eq.® 
and boundary conditions is 

/(0) = | >n =-^(-e x + e 1 ,-e,). (12) 

One can calculate the circulation integral and confirm 
that § dr ■ v s = \ '^ L - A 1/3-quantum vortex is always 
superimposed with a spin defect and is a unique com- 
posite excitation in coherent condensates. When going 
around this vortex, effectively the local triad xyz. is per- 
muted cyclically: xyz — > yzx; this is accompanied by a 
180° rotation around the y axis (see Fig. (|2f ) ) . The sta- 
bility can be studied by examining the core structure. 
Similar discussions were presented for half- vortices [l|| ■ 

When 6l > x and cl < 0, x satisfies trx 2 = 1/2, x = 
X* which corresponds to biaxial nematics [lit ]. A real 
solution is again invariant under the subgroup in Eq.®; 
the submanifold should be isomorphic to M. n = x S 1 , 
Q is the quaternion group introduced before [20(. There 
are quaternion spin defects but no 1/3-quantum vortices 
in this case. The quaternion vortices are analo gou s to 
disclinations in cholesteric liquid crystals 

In conclusion, we have studied discrete symmetries in 
condensates of F = 2 cold atoms and investigated 1/3- 
quantum vortices and pure spin defects. It is worth 
remarking that 1/3-quantum vortices are also natural 
topological excitations in condensates of singlets of three 
atoms (an analogue of Cooper pairs in a two-body chan- 
nel) or condensates of trimers. Rotationally invariant 
Mott states of trimers were pointed out by the authors 
recently [l5j|. If a trimer condensate does exist in optical 
lattices (say created by removing atoms from a trimer 
Mott state), there will then be simple 1/3-quantum vor- 
tices that are featureless in the spin subspace. Note 



that 1 /3-quantum vortices in a cyclic phase on the other 
hand have very rich spin structure as discussed above. 
The one-third circulation quantum can be studied by 
observing fractionalized values in the energy splitting 
of collective modes in experiments similar to those in 
Ref. [H, 0] . FZ is in part supported by the office of the 
Dean of Science at the University of British Columbia, 
NSERC (Canada), Canadian Institute for Advanced Re- 
search and the Alfred P. Sloan foundation. 
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